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ABSTRACT 

We report the detection of hard X-ray emission components in the spectra of six 
nearby, giant eUiptical galaxies observed with the ASCA satellite. The systems studied, 
which exhibit strong dynamical evidence for supermassive black holes in their nuclei, 
are M87, NGC 1399 and NGC 4696 (the dominant galaxies of the Virgo, Fornax 
and Centaurus clusters, respectively) and NGC 4472, 4636 and 4649 (three further 
giant ellipticals in the Virgo cluster). The ASCA data for all six sources provide 
clear evidence for hard emission components, which can be parameterized by power- 
law models with photon indices in the range F = 0.6 — 1.5 (mean value 1.2) and 
intrinsic 1 — 10 keV luminosities of 2 x lO'*" — 2 x 10^^ ergs~^ . Our results imply the 
identification of a new class of accreting X-ray source, with X-ray spectra significantly 
harder than those of binary X-ray sources, Seyfert nuclei or low luminosity AGN, 
and bolometric luminosities relatively dominated by their X-ray emission. We discuss 
various possible origins for the hard X-ray emission and argue that it is most likely 
to be due to accretion onto the central supermassive black holes, via low-radiative 
efficiency accretion coupled with strong outflows. In the case of M87, our detected 
power-law flux is in good agreement with a previously-reported measurement from 
ROSAT High Resolution Imager observations, which were able to resolve the jet from 
the nuclear X-ray emission components. We conflrm previous results showing that the 
use of multiphase models in the analysis of the ASCA data leads to determinations of 
approximately solar emission-weighted metallicities for the X-ray gas in the galaxies. 
We also present results on the individual element abundances in NGC 4636. 

Key words: accretion, accretion disks - galaxies: abundances - galaxies: active - 
galaxies: nuclei - X-rays: galaxies ~ X-rays: ISM 



1 INTRODUCTION 

Observations of nearby galaxies provide convincing evidence 
for the existence of supermassive black holes. Although most 
such galaxies exhibit little or no nuclear activity, dynami- 
cal arguments based on the observed stellar and gas distri- 
butions firmly imply the presence of supermassive compact 
objects in their cores {e.g. Kormendy & Richstone 1995, 
Magorrian et at 1998, Ho 1998, van der Marel 1999). In- 
terestingly, these studies show that virtually all early-type 
galaxies host black holes with masses in the range 10*— a 
few 10^ Mq. In disk galaxies, however, the (rather sparser) 
observations indicate that the black hole masses are typi- 
cally of order 10^ Mq {e.g. SgrA* in our own Galaxy, M31, 
and the maser galaxy NGC 4258), consistent with the black 
hole masses in galaxies being roughly proportional to the 
masses of the spheroidal components. 

Although the central black hole masses in early 
type galaxies are large enough to be the remnants of 
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QSO phenomena (with luminosities ranging from 10'''' to 
lO^^ergs"^ ) and giant ellipticals are known to host radio 
galaxies and radio-loud quasars at high redshifts, in nearby 
ellipticals, only low-luminosity radio cores are commonly ob- 
served (Sadler, Jenkins & Kotanji 1989; Wrobel & Heer- 
shen 1991). Spiral galaxies, in contrast, frequently exhibit 
nuclear activity {e.g. Seyfert nuclei), with optical/UV emis- 
sion, strong, power-law X-ray continua (with a typical intrin- 
sic photon index, F ~ 2) and complex variability observed. 
Such phenomena are usually attributed to standard, thin- 
disk accretion, coupled with a hot coronal plasma with high 
radiative efficiencies. 

Postulating that the early-type systems are simply 
'starved' of fuel to power their activity appears implausi- 
ble since these galaxies possess extensive hot, gaseous halos 
which will inevitably accrete onto their central black holes 
at rates which may be estimated using the Bondi (1952) for- 
mula (a lower limit). This accretion should give rise to far 
more activity than is observed, if the radiative efficiency is 
~ 10 per cent, as is postulated in standard accretion the- 
ory {e.g. Fabian & Canizares 1988). All nearby elliptical 
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galaxies should then host active nuclei with luminosities ex- 
ceeding 10*^ ergs~^ . However, the available data show that 
the total X-ray luminosities of these galaxies rarely exceed 
10*^ ergs^^ , with only a small fraction of this flux being due 
to a central point source. 

Accretion with such high radiative efficiency need not 
be universal, however. As suggested by several authors (Rees 
et al. 1982; Begelman, Blandford & Rees 1984; Fabian & 
Rees 1995) and successfully applied to a number of giant el- 
lipticals in the Virgo cluster (Reynolds et al. 1996; Di Matteo 
& Fabian 1997a; Mahadevan 1997; Di Matteo et al. 1999a), 
the final stages of accretion in elliptical galaxies may occur 
via an advection-dominated accretion flow (ADAF; Narayan 
& Yi 1995, Abramowicz et al. 1995) at roughly the Bondi 
rates. Within the context of such an accretion mode, the 
quiescence of the nuclei in these systems is not surprising; 
when the accretion rate is low, the radiative efficiency of the 
accreting (low density) material will also be low. Other fac- 
tors may also contribute to the low luminosities observed. 
As discussed by Blandford & Begelman (1999; and empha- 
sized observationally by Di Matteo et al. 1999a), winds may 
transport energy, angular momentum and mass out of the 
accretion flows, resulting in only a small fraction of the mate- 
rial supplied at large radii actually accreting onto the central 
black holes. 

If the accretion from the hot interstellar medium in 
elliptical galaxies (which should have relatively low angu- 
lar momentum) proceeds directly into the hot, advection- 
dominated regime, and low-efficiency accretion is coupled 
with outflows (Di Matteo et al. 1999a), the question arises 
of whether any of the material entering into the accretion 
flows at large radii actually reaches the central black holes. 
The present observational data generally provide little or no 
evidence for detectable optical, UV or X-ray emission asso- 
ciated with the nuclear regions of these galaxies. 

In this paper we present a detailed, multiphase analysis 
of high-quality ASCA X-ray spectra for six nearby elliptical 
galaxies; three central cluster galaxies and three giant ellipti- 
cals in the Virgo Cluster. We obtain clear detections of hard 
(mean- weighted photon index, F = 1.2) power- law, emission 
components in the integrated X-ray spectra of all six galax- 
ies. Although our data do not allow us to unambiguously 
identify these components with the nuclear regions of the 
galaxies, we show that this emission is likely to be related 
to the accretion process and the presence of the central, su- 
permassivc black holes, which also power radio activity in 
the galax;ies at varying levels. 

The presence of hard X-ray emission components in the 
ASCA spectra of elliptical galaxies has previously been re- 
ported by a number of authors {e.g. Matsushita et al. 1994; 
Matsumoto et al. 1997; Buote & Fabian 1998). In general, 
however, these studies have associated the observed, hard 
components with a population of binary X-ray sources, the 
spectra of which can be approximated by brcmsstrahlung 
models with typical temperatures of 4 — 7 keV {e.g. Fab- 
biano, Trinchicri & Van Speybroeck 1987; White, Stella & 
Parmar 1988; Makishima et al. 1989; Tanaka et al. 1996; 
Christian & Swank 1998). The spectra of the hard X-ray 
emission components detected in the present sample of el- 
liptical galaxies are significantly harder than those of binary 
X-ray sources or typical AGN, identifying these objects as, 
potentially, a new class of accreting source. (The 'canonical' 



photon index for AGN has F ~ 1.8 - 2; Nandra et al. 1997, 
Reynolds 1997; recent ASCA studies of low-luminosity AGN 
indicate photon indices in the range F ~ 1.6— 1.8; e.g. M81; 
Ishiaki et al. 1996, Turner et al. 1996 or NGC 4258; Mak- 
ishima et al. 1994) The X-ray luminosities of the present 
systems (Lx,i-io = 2 x 10*'^ - 2 x 10"'^ crgs"^ ) are simi- 
lar to those of low luminosity AGN {e.g. Ho, Filippenko & 
Sargent 1997), although the large (~ lO'' Mq) black hole 
masses in their nuclei identifies them as far more inefficient 
radiators. The spectral energy distributions for the galaxies 
in our sample indicate that a significant fraction of their lu- 
minosities are emitted at X-ray wavelengths, with relatively 
low levels of optical and UV emission, which are dominant in 
typical AGN. These data, coupled with the presence of com- 
pact (possibly thermal synchrotron) radio cores, suggest a 
possible, ubiquitous presence of low-level nuclear activity in 
the nearby universe and provide important new constraints 
on the dominant accretion mechanisms in elliptical galaxy 
cores. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target selection 

The targets selected for investigation are six well-studied, 
nearby elliptical galaxies, with high-quality ASCA observa- 
tions available in public archives. The objects include three 
dominant cluster galax;ies; the well known, low-luminosity 
radio galaxy M87 (NGC 4486) at the centre of the Virgo 
Cluster; NGC 1399, the giant elliptical at the centre of 
the Fornax Cluster and NGC 4696, the dominant galaxy of 
the Centaurus Cluster. These three galaxies, and especially 
M87, are key examples of quiescent active nuclei, contain- 
ing nuclear black hole masses of 3 - 5 X 10^ Mo (Ford et 
al. 1995; Harms et al. 1994; Macchetto et al. 1998; Magor- 
rian et al. 1998; although no direct mass measurement for 
NGC 4696 has been made) and exhibiting low-luminosity 
relativistic jets and FR-I-type radio emission. However, the 
luminosities of their nuclei are approximately three orders of 
magnitude less than would be expected if they were accret- 
ing at close to their Bondi rates, with a standard radiative 
efficiency {e.g. Reynolds et al. 1996). 

The sample also includes three other giant ellipticals 
in the Virgo cluster; NGC 4649, 4472 and 4636. These 
galaxies contain almost completely inactive black holes, 
with measured masses ranging from a few 10* M© to a few 
10^ Mq (Magorrian et al. 1998) and predicted X-ray lumi- 
nosities (if the accretion occurs at the Bondi rate with a 
standard radiative efficiency) four to five orders of magni- 
tude greater than is observed (Di Matteo et al. 1999a). 

The central cluster galaxies provide an extreme illus- 
tration of the phenomenon of quiescent black holes. Beside 
possessing FR-I-type radio sources and very large black hole 
masses, these galaxies exist in extremely gas-rich environ- 
ments i.e. in cooling flows at the centres of clusters, and 
are therefore ideal sources in which to study the physics of 
low-radiative efficiency accretion. However, these same rea- 
sons also imply that these galaxies arc by no means typical. 
Despite exhibiting lower luminosities, the properties of the 
other three ellipticals in our sample, which do not exist in 
such preferential locations, may be more easily generalized 
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Table 1. Summajry of the ASCA Observations 



Oluster 






ou 


til 






M87 


18.0 


1993 Jun 07 


13633 


14668 


16874 


16874 


NGC 1399(1) 


29.0 


1993 Jul 15 


16636 


17678 


19646 


19646 


NGC 1399(2) 




1993 Jul 16 


14266 


16961 


19741 


19785 


NGC 4696 


62.6 


1995 May 19 




67650 


70981 




NGC 4472(1) 


18.0 


1993 Jun 30 


17022 


18314 


22242 


22240 


NGC 4472(2) 




1993 Jul 04 


16461 


18781 


22739 


22731 


NGC 4636 


18.0 


1995 Dec 28 


244180 


245450 


183290 


183280 


NGC 4649(1) 


18.0 


1994 Jan 07 


21438 


10669 


40928 


40926 


NGC 1619(2) 


18.0 


1991 .Ian 07 


13775 









Notes: Column 2 lists the luminosity distances (-Dl) to the galaxies in Mpc. For the Virgo Cluster galaxies wc assume a fixed distance 
of 18Mpc. For NGC 1399 and 4696 the distances arc calculated using the optically-determined redshifts and an assumed cosmology of 
-ffo=50 kms^-*^ Mpc^^, Q = 1 and A = 0. Column 3 lists the dates of the observations. Columns 4 — 7 summarize the net exposure times 
(in seconds) in the four ASCA detectors, after all screening and cleaning procedures were carried out. For NGC 4696 only data from 
the SI and G2 detectors were included in the analysis, due to residual calibration errors in the data from the other detectors for that 
observation. For NGC 1399, 4472 and 4649, the observations were made in two parts. For the second NGC 4649 observation, only the SO 
data were of sufficient quality to be of use in the analysis. 



to other systems. These three galaxies have also been stud- 
ied at high radio frequencies and sub-mm wavelengths (Di 
Matteo et al. 1999a) and, when coupled with the reduced 
jet-radio-activity in these systems (and the ASCA data pre- 
sented in this paper), provide further crucial constraints on 
the primary accretion mechanisms in the nuclei of elliptical 
galaxies. 



2.2 The ASCA Observations and data reduction 

The ASCA (Tanaka, Inouo & Holt 1994) observations were 
made over a two-and-a-half year period between 1993 June 
and 1995 December. The ASCA X-ray telescope array 
(XRT) consists of four nested-foil telescopes, each focussed 
onto one of four detectors; two X-ray CCD cameras, the 
Solid-state Imaging Spectrometers (SIS; SO and SI) and 
two Geis scintillation Imaging Spectrometers (GIS; G2 and 
G3). The XRT provides a spatial resolution of ~ 3 ar- 
cmin (Half Power Diameter) in the energy range 0.5 — 10 
keV. The SIS detectors provide excellent spectral resolution 
[AE/E = 0.02(£;/5.9keV)"°-^] over a 22 x 22 arcmin^ field 
of view. The GIS detectors provide poorer energy resolution 
[AE/E = 0.08(£;/5.9keV)"°-^] but cover a larger circular 
field of view of ~ 50 arcmin diameter. 

For our analysis of the ASCA data wc have used the 
screened event lists from the rev2 processing of the data 
sets available on the Goddard-Space Flight Center (GSFC) 
ASCA archive (for details see the ASCA Data Reduction 
Guide, published by GSFC.) The data were reduced using 
the FTOOLS software (version 4.1) from within the XS- 
ELECT environment (version 1.4). Further data-cleaning 
procedures as recommended in the ASCA Data Reduction 
Guide, including appropriate grade selection, gain correc- 
tions and manual screening based on the individual instru- 
ment light curves, were followed. A full summary of the ob- 
servations is given in Table 1. 

Spectra were extracted from all four ASCA detectors. 



except for NGC 4696 for which the SO and G3 data exhib- 
ited residual calibration errors and so were excluded from 
the analysis. The spectra were extracted from circular re- 
gions, centred on the peaks of the X-ray emission from the 
galaxies. For the SIS data, the radii of the regions studied 
were selected to minimize the number of chip boundaries 
crossed (thereby minimizing the systematic uncertainties in- 
troduced by such crossings) whilst covering as large a region 
of the galaxies as possible. Data from the regions between 
the chips were masked out and excluded. The final extrac- 
tion radii for the SIS data are summarized in Table 2. Also 
noted are the chip modes used for the observations and the 
number of chips from which the extracted data were drawn. 
For the GIS data, a fixed extraction radius of 6 arcmin was 
adopted. 

Background subtraction was carried out using the 
'blank sky' observations of high Galactic latitude fields com- 
piled during the performance verification stage of the ASCA 
mission. The background spectra were screened and grade 
selected in the same maimer as the target observations and 
extracted firom the same regions of the detectors. For the sys- 
tems observed in 4-CCD mode, additional SIS background 
spectra were extracted from regions of the chips free from 
significant source counts (we note that this was not possible 
for the M87 observation due to the extended cluster emission 
which covers the full area of the detectors) . The use of these 
on-chip backgrounds in place of the blank-sky backgrounds 
did not significantly alter the results. 

For the SIS data, response matrices were generated us- 
ing the FTOOLS SISRMG software (version 1.1). Where the 
extracted spectra covered more than a single chip, individ- 
ual response matrices were created for each chip, which were 
then combined to form a counts-weighted mean matrix. For 
the GIS analysis, the response matrices issued by GSFC on 
1995 March 6 were used. 
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Table 2. Extraction radii and chip modes for the SIS observations 



Cluster SO SI Chip Mode 

(amin/kpc) (amin/kpc) 



M87 


4.9/25.5 


4.1/21.3 


4(2) 


NGC 1399(1) 


4.2/35.1 


3.7/30.9 


4(1) 


NGC 1399(2) 


4.3/35.9 


4.2/35.1 


4(1) 


NGC 4696 




3.7/66.0 


1(1) 


NGC 4472(1) 


3.7/19.3 


3.7/19.3 


4(1) 


NGC 4472(2) 


3.2/16.7 


3.4/17.7 


4(1) 


NGC 4636 


4.4/22.9 


3.6/18.7 


1(1) 


NGC 4649(1) 


3.8/19.8 


2.9/15.1 


4(2) 


NGC 4649(2) 


3.8/19.8 




2(2) 



The radii of the circular extraction regions for the SIS spectra (in 
aJTcmin and kpc) and the chip modes used in the observations (cither 
1,2 or 4-CCD mode). The numbers in parentheses indicate the num- 
ber of chips contributing to the extracted spectra. For the GIS data, 
a fixed extraction radius of 6 arcmin was used. 

3 THE ANALYSIS OF THE ASCA DATA 
3.1 The spectral models 

The modeling of the X-ray spectra has been carried out using 
the XSPEC spectral fitting package (version 10.0; Arnaud 
1996). For the SIS data, only counts in pulse height analyser 
(PHA) channels corresponding to energies between 0.6 and 
10.0 keV were included in the analysis (the energy range 
over which the calibration of the SIS instruments is best- 
understood). For the GIS data, only counts in the energy 
range 1.0 — 10.0 keV were used. The spectra were grouped 
before fitting to ensure that statistics could be reliably 
used (after background subtraction). 

The spectra have been modeled using the plasma codes 
of Kaastra & Mewe (1993; incorporating the Fe L calcula- 
tions of Liedahl, Osterheld & Goldstein 1995) and the pho- 
toelectric absorption models of Balucinska-Church & Mc- 
Cammon (1992). The data from all four detectors were in- 
cluded in the analysis, with the fit parameters linked to take 
the same values across the data sets. The exceptions to this 
were the emission measures of the (hot) plasma components, 
which model the extended X-ray halos of the galaxies and 
which, due to the different extraction radii used for the dif- 
ferent detectors, were maintained as independent fit param- 
eters. 

The spectra were examined with a scries of spectral 
models. (We have adopted the naming convention of Allen 
et al. 1999 from their analysis of nearby cluster spectra). 
The first model. Model B, consists of an isothermal plasma 
of temperature, kT, and metallicity, Z, in coUisional equi- 
librium, at the optically-determined redshift for the galaxy, 
and absorbed by a column density A'h- Metallicitics are 
measured relative to solar photospheric values of Anders 
& Grevesse (1989) with the various elements assumed to 
be present in their solar ratios. The second model, model D, 
included an additional, cooler plasma component of temper- 
ature, kTi, the normalization of which was a further free pa- 
rameter in the fits. (This cooler component, where present. 



is typically associated with the presence of a cooling flow in 
a cluster or galaxy). The metallicity of the cooler component 
was linked to that of the hotter gas. The cooler component 
was also assumed to be absorbed by an intrinsic column den- 
sity, AA^H, of cold gas, which was a further free parameter 
in the fits. The abundances of metals in the absorbing gas 
were fixed at their solar values (Anders & Grevesse 1989). 

Where it provided a significant statistical improvement, 
a third spectral model, model F, was also examined which 
was similar to model D but with the abundances of various 
individual elements also included as free parameters in the 
fits. This model was statistically preferred over model D for 
M87, NGC 4696 and NGC 4636. 

Finally, a fourth spectral model, Model G, was investi- 
gated in which a power-law emission component was intro- 
duced into the previous best-fitting two-temperature plasma 
model (either model D or F, respectively). The power- law 
emission component was assumed to be absorbed by the 
same intrinsic column density acting on the cooler plasma 
emission component (since both components are expected 
to arise primarily from the central regions of the galaxies 
and/or clusters). 

3.2 The requirement for multi-phase models and 
the metallicity of the X-ray gas 

The results from the spectral modeling are summarized in 
Table 3. We see that the two-temperature model, model D, 
invariably provides a much better description of the ASCA 
spectra for the galaxies than the single-temperature model 
(model B), with a typical reduction in of a few hundred 
for the introduction of only three additional degrees of free- 
dom in the fit. For NGC 1399, 4472 and 4649, the use of 
the two temperature model results in a significant increase 
in the inferred metallicity of the X-ray gas in the galaxies, 
from values of 0.2 — 0.4 solar with spectral model B, to val- 
ues consistent with (or slightly exceeding) the solar value, 
with spectral model D. This is in agreement with the re- 
sults previously-reported by Buote & Fabian (1998; see also 
Buote 1999). 

For M87, NGC 4696 and NGC 4636, we found that 
the fits were further significantly improved by allowing the 
abundances of various individual elements to be included as 
additional free parameters in the fits. (With spectral model 
D, all elements are linked to vary in the same ratio, rela- 
tive to their solar values). Improvements in the measured 
values of several hundred were obtained for these objects 
by allowing the abundances of Mg, Si and S to be included as 
free parameters in the fits (with only a single extra degree of 
freedom being associated with each extra element included 
as a free fit parameter). The results on the individual ele- 
ment abundances for M87 and NGC 4696 are discussed by 
Allen et al. (1999). The results for NGC 4636 are detailed 
in Section 6. 

3.3 The requirement for the power-law 
components 

The results on the power-law components detected in the 
ASCA spectra are summarized in Table 4. In all cases 
the introduction of the power-law component into the two- 
temperature models (models D and F) leads to a highly 
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Table 3. Spectral results for the central cluster galaxies 





Parameters 


Model B 


Model D 


Model F 


Model G 


M87 


kT 

Nh 

kT2 

r 


2 02+"'"^ 
^■"^3-0.02 


9 IQ+0.02 

79+0 04 
"•"'-0.03 

°-"-0.8 


9 17+0.03 
^^■-^'-0.03 
76+0-0-1 

"• '"-0.03 
27+0-07 
"■^'-0.07 
92+0-07 

3-'°-1.09 


2 01+"'"'' 

^•"^-0.03 
"■'^-0.04 

79+0-05 
IJ- '»_0.05 

4-»»_o.86 

1 40+0-37 
^■^"-0.46 

137lif 






2507/1152 


1950/1149 


1558/1145 


1468/1143 



NGC 1399 



kT 


1 09+0.02 
-^-■^^-0.03 


1 -7+0.07 
^- ' ' -0.07 




1 K4+0.07 
^•"^-0.06 


Z 




, 9c:+0.22 






Nh 


0.55l°;il 


0.38«;g 




44+0-28 


kT2 




r, 70+0.04 

"•'^°-0.03 




71+0-04 
'J-'-^-0.03 


ANh 








, C7+0.90 


F 










Ai 








5.8719-73 




12S0/792 


<)ll.S/78<) 




,S7S.(i/787 



NGC 4696 



kT 


9 7^+0.02 
^•'"-0.03 


q ,9+0.10 

•^••5^-0.10 


, 07+0.13 

3-°'-0.21 


, 97+0.26 
•^•^'-0.12 


Z 


1 20+°°^ 
-^■^^-0.04 


97+0-05 


p, 04+0.02 


0.861;;:- 


iVii 


1 09+°°^ 

l.U»_g (J4 


1 44+0-06 
■^•^*-o.ii 


0.89+0 


U.»0_(, 2(i 


kT2 




1 40+0 02 
^•^"-0.04 


1 C9+0.'04 
J^-O^-0.05 


1 46+0-07 


ANh 




00+0-36 

"•""-0.00 


71+0.53 
"•'^-0.60 


59+i °i 

U.O»_g 54 


F 








76+0-59 
"■ '0-0.65 


Ai 








^■'■'J_ig 4 


xVdof 


2124/755 


1219/752 


990.7/748 


954.5/746 



The best-fit parameter values and 90 per cent (Ax^ = 2.71) confidence limits from the analysis of the ASCA spectra for the central cluster 
galaxies. Temperatures (kT), metallicities (Z), column densities (Nh), intrinsic column densities (AATjj), photon indices (F), power- law 
normalizations {Ai), and the normalizations of the cooler plasma emission components in spectral models D, F, G were linked to take the 
same values in all detectors. Only the normalizations of the hotter plasma emission components were allowed to vary independently for 
each detector. Temperatures are quoted in keV and metallicities as a fraction of the solar photospheric -value (Anders & Grevesse 1989). 
Column densities and intrinsic column densities are in units of 10^^ atom cm~^. Power-law normalizations are in units of 10~® photon 
keV-^cm-^s-i at 1 keV. 



significant improvement in the fit. For guidance, a reduction 
in oi Ax^ 10 with the introduction of the power-law 
component (2 extra fit parameters) is significant at approx- 
imately the 99 per cent confidence level (for a fit with 1000 
degrees of freedom and a reduced x^/DOF value ~ 1.0). The 
observed improvements range from Ax^ ~ 30 to Ax^ ~ 600. 



associated with Scyfert nuclei. Figure 1 shows the intrinsic 
1 — 10 keV luminosities as a function of the observed photon 
index. No clear correlation between the two parameters is 
observed. Fig. 2 shows the ASCA spectrum and best-fitting 
model for M87 (spectral model G). For clarity, only the re- 
sults for the SO detector are shown. 



In all cases the slopes of the power-law components are 
significantly flatter (F = 0.6 — 1.5) than the canonical value 
of r ~ 1.8 obtained for Seyfert galaxies {e.g. Nandra et al. 
1997, Reynolds 1997). The weighted mean best- fit photon 
index for the sources in our sample is F = 1.22. The ob- 
served 2-10 keV fiuxes range from 5.7 x 10~" to 8.7 x 10"^^ 
erg cm~^ s~^ , and the intrinsic 1 — 10 keV luminosities of the 
power-law components (corrected for Galactic and intrinsic 
absorption as determined with spectral model G, and calcu- 
lated using the luminosity distances listed in Table 1) range 
from 2.6 x 10*° (NGC 4649) to 2.2 x 10*^ erg s"^ (NGC 
4696). These luminosities are also less than those typically 



The detection of hard, power-hnv omission components 
from all six galaxies is not easily explained as an artifact of 
the analysis. Applying the same method to an ASCA ob- 
servation of NGC 1275, the dominant galaxy of the nearby 
{z = 0.0183) Perseus Cluster, which is known to contain an 
active nucleus, we dotormino a photon index for the nuclear 
emission of F = 2.05 ±0.05, in good agreement with the typ- 
ical values determined for such galaxies (Nandra et al. 1997; 
Reynolds 1997; Turner et al. 1997). The 2-10 koV flux asso- 
ciated with the nucleus of NGC 1275, Fx,2-io ^ 1-8 x 10"^° 
ergcm~^s~^ (implying an intrinsic 1 — 10 keV luminosity 
of ~ 4 X 10** ergs~^ ), is also significantly larger than the 
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Table 3. Spectral results for the Virgo ellipticals 



Parameters 



Model B 



Model D 



Model F 



Model G 



NGC 4472 



kT 




1 66+°-i2 

1-"JD_0.10 




-I qi+0.09 

^•'^^-0.06 


Z 


0.261^2 


1 34+0.33 






Nh 










kT2 




-0.03 




'-^^ ' ^-0.03 






1 10+1" 






r 


















3.82l|-J® 


;2/DOF 


1006/375 


536.7/372 




474.0/370 



NGC 4636 



kT 
Z 
Nh 
kT2 

r 

Ai 

X-/DO]: 



+0.005 
-0.006 
+0.02 
-0.02 

1 42+0-10 



0.661 
0.32 



0.723 
0.35 
1.50 

0.550 
0.00 



+0.030 
-0.021 
+0.03 
-0.02 
+0.12 
0.11 
+0.041 
0.053 
+0.03 
0.00 



g 7Q3 + 0.014 

,+0.03 
'-0.04 
^+0.11 
-0.13 



0.43: 
1.24^ 



0.509 
0.00 



+0.038 
-0.046 
+0.23 
0.00 



0.677 
0.62 
1.03 

0.387 
2.45 



+0.010 

0.009 
+0.13 
0.09 
0.26 
0.25 
+0.052 
0.057 
+2.00 
1.47 
+0.26 



1512/2:55 



llll/2:-!2 



10()0/22<) 



1 4Q 

-^■^^-0.25 

175.0/227 



NGC 4649 



kT 


96+0-0* 

"•^"-0.04 


r, c-,+0.23 
^•^1-0.20 




1 7P;+0.22 
l-'°-0.26 


Z 




^•"^-0.28 




66+0-32 
"•""-0.21 


Nh 


o.8it°:f7 


00+0-13 

U.UU_() QQ 




00+0-26 
"•""-0.00 


kT2 




79+0 04 

"■"'-0.03 




n 77+0.05 
"•' '-0.04 


ANh 




9 1 1+0.62 
^•-^^-0.70 




0/1+0.66 
^••3"*_0.80 


r 








"•^^^-0.54 










9 94+5 52 
^•''^-1.45 


xVdof 


606.3/240 


318.8/237 




291.1/235 



fluxes associated with the harder power-law components de- 
tected in the present sample of objects. Secondly, when we 
apply the same analysis method (using spectral model D) 
to an observation of the central regions of the Coma Clus- 
ter (2 = 0.0232), where the observations are centred on 
the X-ray centroid of the cluster rather than an individ- 
ual galaxy (the cluster does not contain a single, optically- 
dominant galaxy in its core), wc find no improvement in the 
fit (Ax^ = 0.0) on the introduction of a hard, power-law 
component. The upper (90 per cent confidence) limit to the 
2 — 10 keV flux of any power-law component with a photon 
index T = 2.0 in the ASCA data for centre of the Coma 
Cluster using spectral model D is Fx. 2-10 < 2.2 x 10^^^ 
ergcm~2s-i^ significantly less than the measured values 
for M87 and NGC 4696. We also note that the applica- 
tion of a similar analysis procedure to ASCA observations of 
more distant, luminous cooling-flow clusters docs not, gen- 
erally, indicate the presence of hard (F ~ 1.2) power-law 
components in these systems. (We do not expect to detect 
hard, power-law components with luminosities comparable 
to those reported here in clusters an order of magnitude 
(or more) more X-ray luminous than the Virgo or Centau- 
rus clusters; c.f. Section 5.3) This further suggests that the 
hard, power-law components detected in the present sample 
of elliptical galaxies are not artifacts due to systematic errors 



in the modeling of difi^use cluster/galaxy emission in these 
systems. (A more detailed discussion of the ASCA data for 
the Perseus and Coma clusters is presented by Allen et al. 
1999. For a discussion of more distant, luminous cooling-flow 
clusters see Allen 1999). 

Finally, we note that the data for NGC 4636 exhibit a 
systematic discrepancy at low energies (E < IkeV) between 
the SO and SI detectors, which contributes signiflcantly to 
the measured . (This discrepancy is also noted by Buote 
1999). Ignoring the data in this region and repeating the 
analysis with spectral model G leads to a lower x2 value 
(x^ = 284 for 201 degrees of freedom) and constraints on the 
power-law component in good agreement with those listed 
in Table 2. 



3.4 A bremsstrahlung model for the hard X-ray 
emission 

We have investigated whether the hard X-ray emission from 
the elliptical galaxies can also be parameterized by a sim- 
ple bremsstrahlung model. To do this, the power-law com- 
ponent in spectral model G was replaced with a thermal 
bremsstrahlung component, with the temperature, fcTbrem, 
and normalization, j42, as flt parameters. The results from 
the fits with this modified model are summarized in Table 
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Figure 1. The intrinsic 1 — 10 keV luminosities of the detected 

power-law components as a function of the observed photon index, 
r. Error bars are 90 per cent (Ax^ = 2.71) confidence limits. No 
correlation between the luminosity and photon index is observed. 

5. The values obtained indicate that the bremsstrahlung 
model provides as good a parameterization for the hard 
X-ray emission from the galaxies as the power-law model. 
Generally, the results constrain the temperatures of the 
bremsstrahlung components to be > 10 keV, with the upper 
limits on the temperatures unconstrained. (ASCA data in 
the 0.6—10.0 keV band cannot reliably constrain the temper- 
atures of plasmas hotter than ^ 10 kcV). The results on the 
other fit parameters are essentially identical to those listed in 
Table 3. The fluxes and intrinsic luminosities associated with 
the bremsstrahlung components are also in good agreement 
with the values listed in Table 4, using the power-law model. 
The measured temperatures for the bremsstrahlung compo- 
nents are significantly hotter than those of typical Galactic 
binary sources [e.g. Christian & Swank 1998) and support 
the interpretation for the origin of the hard X-ray emission 
from the galaxies in terms of low radiative efficiency accre- 
tion flows, given in Section 4. 

We have examined whether the statistical quality of the 
fits with spectral model G (or the modified model G, incor- 
porating the bremsstrahlung component) can be improved 
by the addition of further model components. In all cases, 
the introduction of a further power-law or bremsstrahlung 
component did not significantly improve the fits. The intro- 
duction of a third, thermal plasma component (as described 
in Section 3.1) also leads to no significant improvement (at 
the 99 per cent confidence level) for any object, except NGC 
4696. For this source, the introduction of a third plasma 
component (with the element abundances linked to be equal 
to those in the other two components) leads to a drop in 
of Ax^ = 37, and slightly modified constraints on the power- 
law emission (F = 1.10+o;44 and Ai = 5.8tl''t x 10^"* photon 
keV~^cm^^s^^). We note that the relatively high temper- 
ature {kT ~ 3.3 keV) for the hotter thermal component in 



NGC 4696 simply reflects the virial temperature of the host 
cluster. NGC 4696 is the most distant galaxy included in 
the present study and the 3.7 arcmin (66 kpc) SI aperture 
used includes significant flux from the extended cluster gas. 

Finally, we note that several of the x^ values listed in 
Tables 3 — 5 (in particular, those for the galaxies with the 
highest signal-to-noise ratios in their ASCA spectra) indi- 
cate that the best- fit models are, formally, statistically un- 
acceptable. However, the high values obtained are pri- 
marily due to residual systematic errors in the instrument 
response matrices and plasma emission models (although we 
do not expect these to significantly effect the main conclu- 
sions reported here.) 



4.1 Accretion processes 

The presence of hard, power-law X-ray emission from astro- 
physical sources provides a discriminating signature of accre- 
tion processes around black holes. Such emission is usually 
attributed to the presence of a hot, tenuous coronal plasma 
(probably magnetically) coupled to an accretion disk or a 
hot (ADAF-type) accretion flow, but can also arise in the 
shock sites of a jet. Other astrophysical processes tend to 
give rise to softer X-ray emission, with a steeper photon in- 
dex (when this emission is modeled as a simple power-law) . 
The flat slopes of the power-law components detected in 
our sample of elliptical galaxies, together with the dynam- 
ical evidence for supermassive black holes with masses of 
10* — 10^° M0 in the nuclei of these objects, argues strongly 
for this emission being intimately related to the accretion 
process. 

The relatively broad point spread function of the ASCA 
mirrors (Section 2.2) does not allow us to resolve the X- 
ray emission from the galaxy cores into more than a single 
integrated spectrum. We therefore cannot separate the X- 
ray emission associated with the jets in these objects {e.g. 
knot A in M87, which previous ROSAT High Resolution Im- 
ager observations have resolved; Reynolds et al. 1996, Harris, 
Biretta & Junor 1997, 1998) from the accretion disks them- 
selves. However, the fact that the galaxy with the strongest 
radio emission, M87, in which the radio and X-ray jet emis- 
sion is relativistically beamed towards us, exhibits a signifi- 
cantly weaker hard X-ray component than NGC 4696 (which 
has a twin lobe radio structure and a total 4.85 GHz radio 
flux ~ 40 times lower than M87) at least suggests that the 
jet emission does not dominate the detected hard X-ray flux. 
This is further supported by the detections of hard, power- 
law X-ray components, with similar characteristic slopes, in 
objects such as NGC 4636, 4649, and 4472, in which the 
radio activity is at the level of a mjy or less. Sharp cut- 
offs observed in the radio spectra of these objects (Di Mat- 
teo et al. 1999a) also strongly constrain the populations of 
non-thermal particles responsible for synchrotron radiation 
and synchrotron self-Compton emission (in X-rays) from a 
jet and/or outflow associated with a low radiative-efficiency 
accretion flow (although contributions from a non-thermal 
distribution of relativistic electrons might still be important 
in the central cluster galaxies with more dominant radio ac- 
tivity; Di Matteo et al. 1999b). 



NGC 4472 
NGC 4649 



THE ORIGIN OF THE HARD X-RAY 
EMISSION 
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4-1.1 The effects of intrinsic absorption 

The ASCA spectra alone cannot reliably discriminate be- 
tween whether the observed power-law components are in- 
trinsically flat or have steeper photon indices which have 
been modified by the effects of intrinsic absorption over and 
above that accounted for by spectral model G. (This un- 
certainty is primarily due to the complex spectrum of the 
diffuse, thermal emission at soft X-ray energies). In general, 
the inclusion of an extra ~ 5 x 10^^ — 10^^ atom cm~^ of 
intrinsic absorption associated with the power-law compo- 
nents leads to power-law photon indices of F ~ 2, and pro- 
vides similarly good fits to the ASCA data (although the 
determination of confidence limits on the fit parameters be- 
comes difficult with such models). Conversely, if the true 
intrinsic column densities acting on the hard X-ray compo- 
nents are overestimated using spectral model G, then the 
intrinsic photon indices of these components will actually 
be flatter than the values listed in Table 4. 

An intrinsic photon index of F ~ 1.8 — 2 is typically 
observed in AGN [e.g. Nandra et al. 1997; Reynolds 1997; 
Turner et al. 1997), where the accretion is thought to occur 
via a thin disk, and where the radiative efficiency is ~ 10 
per cent. However, as discussed in Section 1, if the accretion 
in the present sample of elliptical galaxies were to proceed 
in this manner, the X-ray luminosities associated with their 
power-law emission components should be 3 — 5 orders of 
magnitude larger than the observed values. There is no sim- 
ple way in which to modify the observed X-ray fluxes (par- 
ticularly when one takes into account the joint ASCA and 
ROSAT results; Section 5.1) and measured photon indices 
to agree with the predictions from the standard thin-disk 
accretion models by simply including extra absorption on 
the power-law components. 

IRAS measurements of the 60 and IQO/im fluxes from 
the galaxies (determined using the IPAC SCANPI software 
applied to co-added IRAS scans) also constrain the X-ray 
luminosities that may be absorbed and reprocessed to in- 
frared wavelengths in these systems. For M87, the observed 
60 and lOO/im luminosities (Leo/jm ~ 8±1 x 10*^ ergs~^ and 
iioOfim ~ 4± 1 X lO'*^ ergs"'^ , respectively) are comparable 
to the 1 — 10 keV luminosity of the hard X-ray component 
(Table 4). The infrared luminosities for the other galaxies 
are summarized in Table 6. In general, the observed 100/im 
luminosities do not exceed the 1 — 10 keV values by more 
than an order of magnitude (with the exception of NGC 
4649, for which the observed lOO/im luminosity is ~ 500 
times larger than the 1 — lOkeV value) and the presence of 
an absorbed active nucleus with an X-ray luminosity 4—5 or- 
ders of magnitude more luminous in X-rays than the values 
listed in Table 4 can be firmly ruled out. We note, however, 
that although the X-ray and infrared data argue strongly 
against the observed hard, power-law components being due 
to intrinsically absorbed, steep (F ~ 2) X-ray spectra with 
luminosities ^ 10''® ergs"' , the very flat photon indices ob- 
served in some galaxies (with F < 1) suggest that intrinsic 
absorption may play some role. Such flat slopes could result 
from absorption of few xlO^'atom cm"'^ acting on emis- 
sion spectra with intrinsic photon indices, F ~ 1.4 (which 
could be produced by low-radiative efficiency accretion flows 
with electron temperatures < lOOkeV; Section 4.1.2). This is 
consistent with the observation that M87, which has an X- 



ray photon index, F = 1.4, exhibits a lower LiooMm/ix,i-io 
ratio than those galaxies with flatter X-ray slopes. 



4-1.2 Low radiative efficiency accretion flows 

Although the observed properties of the elliptical nuclei are 
not easily explained within the context of unifled models 
for Seyfert galaxies, the characteristic hard X-ray, radio and 
(lack of) optical emission from these sources (i.e. the typ- 
ically non-thermal character of their spectra) can be ac- 
counted for by models of low-radiative efficiency accretion 
(such as an ADAF; Narayan & Yi 1994, 1995) coupled with 
outflows/winds (Blandford & Begelman 1999), with the ac- 
cretion occurring at approximately the Bondi rates. 

In an advection dominated flow, which occurs at accre- 
tion rates M < a^A/Edd (where M^dd is the accretion rate 
at the Eddington limit) the viscosity (parameterized by the 
constant a > 0.1) is assumed to dissipate most of the energy 
locally into ions. The ions cannot cool (Coulomb scattering 
between the ions and electrons is very inefficient in the low 
density plasma; no other form of electron-ion coupling is as- 
sumed and the gas is two-temperature) and flow inwards, 
carrying an increasing amount of thermal energy. In the ab- 
sence of an energy sink, which would normally be present in 
the form of efficient radiation in a thin disk, the gas will be 
unbound (i.e. have a Bernoulli parameter greater than zero) 
unless the binding energy is transported radially outward 
by the viscous torque, in the form of a wind (Blandford & 
Begelman 1998). Within the context of low-radiative effi- 
ciency accretion fiows (unavoidably coupled with outflows), 
most of the accreted mass, angular momentum and energy 
will be lost at large distances in the flows and the central 
pressures and densities will be much reduced. 

In these more generalized models for low-radiative ef- 
ficiency accretion flows, the density varies as p oc r~^^'^^^ 
(where the accretion rate M cx r^, with p ~ 1) such that 
the emission in the central regions is highly suppressed. The 
thermal electrons, at a temperature of ~ 100 keV, radiate by 
synchrotron emission, inverse Compton scattering (of syn- 
chrotron photons) and bremsstrahlung emission. Most of the 
synchrotron emission and its Comptonization will occur in 
the innermost regions (within a few Schwarschild radii) of 
the flow (where the temperature is the high enough, > 10^ K, 
for these processes to become important) and will therefore 
be highly suppressed in the presence of winds ^ (as em- 
phasized by the observations and modeling discussed by Di 
Matteo et at 1999a,b; see also Quataert & Narayan 1999). 
Bremsstrahlung emission, in contrast, arises from all radii 
in the flow and should, therefore, dominate the emission 

t Although the luminosities due to both bremsstrahlung and syn- 
chrotron processes vary as the square of the accretion rate (i.e. 
the density), in the presence of a wind the accretion rate will be 
much reduced in the inner regions where most of the synchrotron 
emission is produced. More importantly, in the very low density 
inner regions the electron temperature profile is much flatter; the 
synchrotron power oc is highly suppressed (Di Matteo et al 
1999b; Quataert &; Narayan 1999). Comptonization of this com- 
ponent is also highly suppressed, its importance also depending 
on the scattering optical depth r oc M, which always remains low 
in the presence of a wind. 
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Figure 2. (Left) The ASCA SO data and best-fitting two-temperature-plus-power-law model for M87 (with the abundances of various 

individual elements also included as free fit parameters; Allen et al. 1999) and the residuals to the fit in units of x- (Right) The best-fit 
model for the SO data, showing the relative contributions of the hot and cool plasma emission (dotted curves) and power-law components 
(lower solid curve) to the total flux (bold solid curve). 



from sources in which low radiative efficiency accretion is 
associated with winds. The X-ray spectrum due to thermal 

bremsstrahlung in the optically thin gas has a very flat spec- 
trum up to its cut off frequency at about hv ~ kT. Most 
of the bremsstrahlung luminosity, oc p^T~ 



with p oc M^^Mr~^^^^'^ , originates at larger radii where 
the density is the highest. The temperature profile in an 
ADAF with winds is virtually constant with radius with 
kT w 100 - 200 keV for r < a few hundred Schwarschild 
radii, and tracks the virial temperature in the outer re- 
gions. Even in the presence of a very strong wind the 
bremsstrahlung spectrum would not be highly suppressed 
up to energies kT < 10 koV, as the minimum temperature 
in the outer region of the flow is ~ 10^^/routK and the 
outer radius Tout is of the order of a 10^ — 10* Schwarschild 
radii. As discussed by Di Mattco et al. (1999b), the rates at 
which material is required to be fed from the hot interstel- 
lar medium to the outer regions of the accretion flows, to 
explain the luminosities of the observed hard. X-ray com- 
ponents, are consistent with the expected Bondi accretion 
estimates of 0.1 — 1 Moyr^"*^. The omission from such a flow 
would be dominated by a thermal bremsstrahlung at tem- 
peratures < 100 keV (as expected in the outer regions of 
the flows) resembling a hard (F ~ 1.4) power law in the 
(1 - 10 keV) ASCA band (as required by the data). The dif- 
ferences in luminosity between the objects in our sample can 
then be ascribed to differences in the black hole masses and 
Bondi accretion rates, with the hard, power-law components 
in the central cluster galaxies being more luminous due to 
their higher density environments. This is illustrated in Fig. 
3 where we show the intrinsic (1 — 10 keV) luminosities of 
the hard, power-law components as a function of the bolo- 
metric luminosity of the X-ray gas within a radius of 10 kpc 



i/V^exp-'"'/'^^ 



in the galaxies. A clear correlation between these quantities 
is observed, with NGC 4696, the dominant galaxy of the 
Centaurus Cluster and the galaxy with the most luminous 
hard, X-ray component, also exhibiting the largest X-ray gas 
luminosity within a radius of 10 kpc. Detailed modeling and 
discussion of these issues, and of the effects of non-thermal 
particle distributions in the winds and/or jets, are presented 
by Di Matteo et al. (1999b). 



4.2 Other sources of X-ray emission 

Although the observed properties of the hard X-ray compo- 
nents can be accounted for by models of low-radiative effi- 
ciency accretion onto the central supermassivc black holes 
in the galaxies, a variety of other sources (given the large 
field of view of the ASCA instruments) may contribute to 
the integrated X-ray spectra. 



4.2.1 Binary X-ray sources 

Undoubtedly, some contribution towards the harder X-ray 
emission from the galaxies will be made by binary X-ray 
sources {e.g. Canizares, Fabbiano & Trinchieri 1987; Fab- 
biano 1989). Bright Galactic X-ray binaries and black-hole 
candidates exhibit persistent X-ray luminosities in the range 
1036 „ iqSs ergs"^ {e.g. White et al. 1988). The lunfinosi- 
ties associated with the power-law components detected in 
the Virgo ellipticals could therefore be accounted for by a 
few 10^ — 10* such sources (or ~ 10* — 10® sources for the 
more luminous central cluster galaxies). The luminosities of 
the power-law components in NGC 1399, 4472, 4636 and 
4649 are consistent with an extrapolation of the Lx/is re- 
lation determined for irregular and spiral galaxies (Fabbiano 
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& Trinchieri 1985), in which the X-ray emission is dominated 
by discrete sources (binary X-ray sources, stars and central 
active nuclei). Thus, binary X-ray sources could contribute 
at a significant level (comparable to the accretion fiows; see 
below) to the harder X-ray emission from these galaxies. 
For M87 and NGC 4696, however, the X-ray luminosities 
of the hard, power-law components are am order or mag- 
nitude larger than the values estimated from the Lx/Lb 
relation, based on their observed blue luminosities. Binary 
X-ray sources arc therefore unlikely to contribute signifi- 
cantly to the hard, power-law components detected in M87 
and NGC 4696. (The result for M87 is strongly supported 
by the ROSAT HRI observations discussed in Section 5.1). 

In addition, although individual binary sources can ex- 
hibit very hard X-ray spectra, a more typical spectrum in 
the 2 — 10 keV band would have a photon index (using a sim- 
ple power-law model) in the range F = 1.5 — 2.5 or, using a 
thermal bremsstrahlung model, a temperature, kT ~ 4 — 7 
keV {e.g. Fabbiano et al. 1987; White et al. 1988; Mak- 
ishima et al. 1989; Tanaka et al. 1996; Christian & Swank 
1997), with little or no associated intrinsic absorption (Fab- 
biano et al. 1987). These power-law slopes are significantly 
steeper (or, cquivalently, the bremsstrahlung temperatures 
are lower) than the observed values for the hard emission 
components detected in the present sample of galaxies. (We 
recall that the observed photon indices show no clear corre- 
lation with the luminosities of the hard, power-law compo- 
nents; Fig. 1). We have examined the effect of including a 
second power-law component, with a fixed photon index of 
r = 2.0, in the analysis of the ASCA data for NGC 4472 and 
4636 with spectral model G. (We assume that the intrinsic 
absorption acting on both power-law components and the 
cooler plasma component is the same since the ASCA spec- 
tra cannot easily constrain more detailed models.) For both 
systems, we find that the fits are not significantly improved 
by the introduction of the steeper power-law components 
(Ax^ = 0.0; sec also Section 3.4), although the maximum 
(90 per cent confidence) 2 — 10 kcV fluxes associated with 
these components (2.0x10"^^ and 3.3x10"" ergcm"^s"^ ) 
are approximately 30 and 60 per cent of the values listed in 
Table 4, respectively. With the steeper components included 
at their maximum allowed levels, the photon indices of the 
harder power-law components in NGC 4472 and 4636 are 
reduced to 0.71° ^ and 0.9lg g, respectively. 

We note that the relatively low source fiuxes and high 
background count rates at harder energies, together with the 
broad point spread function of the ASCA mirrors, prevent 
us from placing useful constraints, from the ASCA imaging 
data, on the relative contributions to the hard X-ray fluxes 
from central point sources and extended emission compo- 
nents in the galaxies. (For M87 and NGC 4696, the extended 
cluster emission dominates the X-ray spectra across essen- 
tially the entire ASCA energy band; Section 5.3). Future 
observations with the Chandra Observatory will provide an 
answer to this question. 

We conclude that binary X-ray sources are unlikely to 
dominate the hard, power-law emission detected in M87 and 
NGC 4696, although they will contribute to the measured 
fluxes in the galaxies with weaker hard. X-ray components. 
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Figure 3. The intrinsic 1 — 10 keV luminosities of the hard, 

power-law components as a function of the bolomctric luminosity 
of the X-ray gas within a radius of 10 kpc of the galaxy centres. 
(All luminosities are in units of 10^'^ ergs^^ .) A clear correla- 
tion between the quantities is observed, with a best-fitting slope 
(determined from a fit with a power-law model) of between 1 and 
2. 



4.-2.2 Cosmic XRB fluctuations 

It is important to consider the effects of fluctuations in the 
Cosmic XRB on the results, particularly given the close 
match between the weighted mean slope of the detected 
power-law components (F ~ 1.2) and the XRB (F = 1.4). 
Following Barcons, Fabian & Carrera (1998), we have esti- 
mated the confusion noise, Uc (the flux equivalent to a 1 — u 
variation in the XRB intensity histogram) for the ASCA ob- 
servations. Assuming Gc = 2.0 x 10~^^Q,^^^ , where Q. is the 
beam area of the ASCA SIS spectra in degree^ , and normal- 
izing to the GINGA observations of Butcher et al. (1997), 
we obtain the confusion limits listed in column 6 of Table 4. 
For the strongest sources (M87 and NGC 4696) the detected 
fluxes are much larger (a factor 40 — 60) than the confusion 
limits. For even the weakest sources {e.g. NGC 4636), the 
detected fluxes are 4 — 5 times larger than the confusion 
limits. We also recall that the use of an independent on-chip 
background-subtraction method in the analysis of the ASCA 
observations made in 4-CCD mode (section 2) provides very 
similar results on the power-law components. 



4-2.3 Cosmic ray electrons and inverse Compton emission 

A third mechanism that could contribute to the detected 
X-ray fluxes is inverse Compton emission due to primary 
cosmic ray electrons in the intracluster medium. Detailed 
models {e.g. Sarazin 1999) show that for steady particle 
injection (with a given power-law distribution) the inverse 
Compton spectra relax into a steady-state form. Although 
these models predict ubiquitous EUV and soft X-ray emis- 
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Table 4. Properties of the hard, power-law components 



FIT PARAMETERS 
r Ai 



■^-X.l-lO -Pk,2-10 



CON. 
LIMIT 

Cc,2-10 



STATISTICS 
Ax^ xVDOF 



M87 

NGC 1399 
NGC 4696 
NGC 4472 
NGC 4636 
NGC 4649 



1.40 
0.86 
0.76 
0.78 
1.49 
0.55 



+0.37 
0.46 
+0.,57 
0.63 
().,")!) 
0.6."; 
+0.55 
-0.65 
+0.25 
0.25 
+0.72 
0.54 



159 
55 
+9.73 
4.05 
+62.1 
16.4 
+5.88 
-2.57 

2.241S-52 



137 
5.87 
22.0 
3.82 



)+13.7 
-10.6 

10 4+^'^ 
214+99 

3 01+"-2^ 
•^■'^^-o.io 

^+0.14 
-0.10 
a+0.38 
-0.24 



42.2^ 



2.8O: 
2.58"' 



86.9 
9.26 
42.0 
7.09 
5.72 
6.21 



17.0 
15.6 
0.30 
0.23 
+ 12.9 

9.9 
+0.26 
0.08 
+0.16 
0.16 
+0.24 
0.36 



1.52 

1.24 
1.12 
1.05 
1.32 
1.08 



90.0 

63.2 
38.8 
62.7 
585 
27.7 



1468.1/1143 

878.6/787 
954.7/746 
474.0/370 
475.0/227 
291.1/235 



A summary of the results on the hard, power-law emission. Column 2 lists the measured photon indices (F) . The normalizations of the 
power-law components (Ai) at an energy of IkeV are quoted in units of 10~^ photon keV~^cm~^s~^. Luminosities (ix) in the 1 — 10 
keV band are corrected for absorption and quoted in units of lO^^'ergs"^. The observed fluxes in the 2-10 keV band (Fx) are not 
corrected for absorption and are quoted in units of 10~^^ ergcm~^s~^ . <Tc values are the confusion limits for the SO detectors in the 
2 — 10 keV band in units of 10~^^ ergcm~^s~^ . The Ax^ values are the improvements in x^ obtained with the introduction of the 
power-law component into the two-temperature models (either D or F as appropriate). The total x^ and degrees of freedom (DOF) for 
the fits arc also listed. The errors on F and Ai are the 90 per cent confidence limits (Ax^ = 2.71) on a single interesting parameter. The 
errors on the Lx and Fx values account for the joint errors on the photon indices and normalizations. 



Table 5. Additional constrains on the hard X-ray components 
using the bremsstrahlung model 



rem A2 x^/DOF 



M87 


> 9.8 


234^2 


1467.1/1143 


NGC 1399 


> 18.0 


33.0+g7 


879.5/787 


NGC 4696 


> 29.0 


1661^5 


956.0/746 


NGC 4472 


> 18.6 


23.711^ 


475.5/370 


NGC 4636 


> 8.4 


14.6+1-? 


473.1/227 


NGC 4649 


> 26.0 


2i.2it; 


293.3/235 



A summary of the results on the hard X-ray components using the 
bremsstrahlung model. Column 2 lists the measured electron tem- 
peratures (fcThrem) in kcV. Column 3 lists the normalizations {A2) 
in units of 10^'' photou keV~^cm~-^s~^. Column 4 summarizes the 
X^ values (which are similar to those obtained with the power-law 
model in Table 4) and the numbers of degrees of freedom in the fits 
(DOF). Errors are 90 per cent confidence limits (Ax^ = 2.71) on a 
single interesting parameter. 

sion in clusters, due to electrons with 7 --^ .'300 (which have 
the longest loss times and could be reasonably be injected 
since z < 1), at harder energies their contribution is negli- 
gible. Significant hard X-ray emission (1 — 50kcV) can only 
be produced by electrons with 7 ~ 10'^ — 10* i.e. by parti- 
cles with rather short lifetimes (as set by inverse Compton 
losses; Uoss ~ 10^ yr) and would only be present in clusters 
in which substantial particle injection has occurred since 
z < 0.1. Assuming that particles arc accelerated in shocks 
in the intracluster medium, one would only expect diffuse 
hard X-ray emission in clusters undergoing, or having re- 
cently experienced, a major merger event. The Virgo cluster 



docs not show evidence for recent merger activity in its cen- 
tral regions (Allen et al. 1999). In addition, oven in cases 
were such particle injection does occur, the expected steady 
state photon index where inverse Compton losses dominate 
would be r « 2.1, significantly steeper than the observed 
slopes of the power-law conipononts. 

Inverse Compton emission from the radio lobes is a fur- 
ther possible source of power-law X-ray emission from the 
galaxies. Such emission is normally very weak but has been 
detected from the hot-spots of Cygnus-A (Harris, Carilli & 
Perley 1994; Reynolds & Fabian 1996) and Fornax-A (Feigel- 
son et al. 1995; Kaneda et al. 1995). The photon index of the 
resulting X-ray emission should have a photon index similar 
to, or steeper than, that of the radio emission which, for the 
radio lobes, will be significantly steeper than P 1.2. Thus, 
although inverse Compton emission is expected at some level 
in all of the systems, it is unlikely to contribute significantly 
to the hard, X-ray components reported here. 



4.2. 4 X-ray 'reflection' from cold material 

A final possibility for obtaining a relatively flat photon index 

in the ASCA band is the situation where the observed X-ray 
flux is dominated by photons Compton scattered off cold, 
optically-thick material close to the central X-ray source 
[e.g. Lightman & White 1988; George & Fabian 1991; Matt, 
Perola & Piro 1991. This situation also requires that the 
primary X-ray source is heavily obscured, a possibility ar- 
gued against in Section 4.1). The scattered X-ray spectrum 
will include absorption and emission features due to var- 
ious elements in the scattering medium and in particular 
should exhibit a strong, fluorescent Fe-K emission line at 
6.4 koV. ASCA observations of the Circinus galaxy (Matt 
et al. 1996) and NGC 6552 (Reynolds et al. 1994) exhibit 
such 'reflection-dominated' spectra with strong (redshifted) 
6.4 keV emission lines with equivalent widths of ~ 2 keV. 
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Table 6. IRAS luminosities and 6.4 keV line limits 



^lOOjjm i'100;Um/i'X,l- 



E.We., 



M87 
NGC 1399 
NGC 4696 
NGC 4472 
NGC 4636 
NGC 4649 



78 ±9 
< 87 

300 ± 90 
35 ± 12 
31 ±7 

340 ± 10 



40± 10 
81 ± 26 

900 ± 260 
21 ± 17 
< 60 

1260 ± 20 



0.95 ± 0.36 
7.8 ±2.7 
4.2 ±2.1 
7.0 ±5.7 

< 22 
490 ± 60 



< 90eV 

< 150eV 

< 190cV 

< 300eV 

< 130eV 

< 210eV 



A summary of the 60 and lOOfirn IRAS luiiiiiiositics for the galaxies (in units of 10^ 



) determined with the SCANPI software 



and co-added IRAS scans. Error bars are the r.m.s. deviations of the residuals after baseline subtraction. Where no detection of a source 
was made, an upper limit equivalent to 3 times the r.m.s. deviation is quoted. Column 4 lists the ratio of the 100/im luminosities and 
1 — lOkeV luminosities of the hard X-ray components. Column 5 lists the 90 per cent confidence upper limits on the equivalent widths, 
E.W. (in eV, relative to the power-law continua) of an intrinsic, narrow (instrumental resolution) 6.4 keV (Fe K) emission line in the 
ASCA spectra. 



We have searched for the presence of intrinsic 6.4keV emis- 
sion linos in the ASCA spectra for the galaxies. In all cases 
we find no improvement to the fits obtained with spectral 
model G on the introduction of a narrow 6.4 keV emission 
Une (Ax^ = 0.0), and are able to place 90 per cent confidence 
limits {Ax^ = 2.71) on the maximum equivalent widths (rel- 
ative to the power-law continua) of any such lines as listed in 
Table 6. We conclude that the flat, power-law components 
detected in the elliptical galaxies are unlikely to be due to 
this process. 

4.3 Sumrriciry of the results on the origin of the 
hard X-ray components 

We have argued that the hard (F ~ 1.2) power-law, X- 
ray components detected in the ASCA spectra are likely 
to be due to bremsstrahlung emission from low-radiative ef- 
ficiency accretion flows onto the central supermassive black 
holes, coupled with winds/outflows. This is a very different 
situation from Seyfert galaxies, where the power-law X-ray 
emission is normally attributed to thermal inverse Compton 
scattering of the soft, disk radiation field. 

We have examined a range of other mechanisms that 
could contribute towards the detected X-ray fluxes. Some 
contribution towards the harder X-ray emission could be 
made by the jets in these sources, although the absence of 
a clear correlation between the 5 GHz radio and hard X-ray 
luminosities suggests that the jet emission does not domi- 
nate the hard, power-law flux. The close agreement between 
the ASCA and ROSAT (Section 5.1) X-ray fluxes for M87, 
and the IRAS infrared measurements for the galaxies, argue 
strongly against the observed flat photon indices and low 
luminosities being simply due to intrinsic absorption act- 
ing on canonical Seyfert-like spectra (due to accretion at 
the Bondi rates with a standard radiative efficiency). Some 
level of intrinsic absorption is possible, however, and could 
explain the very flat power-law components (F < 1) ob- 
served in some objects (further flattening the spectra of the 
low-radiative efficiency accretion flows, which should have an 
intrinsic photon index in the ASGA band of F ~ 1.4). Limits 
on the 6.4 keV emission-line fluxes rule out a significant con- 



tribution to the detected hard, power-law components fi-om 

X-rays Compton scattered off cold, optically-thick matter 
surrounding the central X-ray sources. 

Binary X-ray sources are likely to contribute to 
the harder X-ray emission from the galaxies with lower- 
luminosity power-law components (the Virgo ellipticals and 
NGC 1399). For M87 and NGC 4696, however, the measured 
power-law luminosities are an order of magnitude greater 
than the values predicted from the Lx/Lb relation for spiral 
and irregular galaxies. Moreover, the typical X-ray spectra 
of X-ray binaries and black hole candidates are significantly 
steeper than the observed hard, power-law components. In- 
verse Compton emission from the radio lobes in the galaxies 
and primary cosmic-ray electrons in the intracluster medium 
should not contribute significantly to the detected power-law 
fluxes and, where present, will typically produce a photon 
index in the ASCA band significantly steeper than the mea- 
sured values. Fluctuations in the Cosmic XRB should not 
significantly effect our results. 

Future X-ray observations at high spatial resolution 
made with the Chandra Observatory will be crucial in estab- 
lishing the contributions from the various emission mecha- 
nisms outlined above and unambiguously identifying the ori- 
gin of the hard X-ray components. Deep X-ray spectroscopy 
with XMM and ASTRO-E will allow us to examine variabil- 
ity in the X-ray emission and search for broad, iron emission 
features associated with the power-law components. If the 
hard X-ray emission indeed originates from bremsstrahlung 
processes in the outer regions of low radiative-efficiency ac- 
cretion flows, then the variability timescales observed should 
be longer than those Eissociated with typical Seyfert nuclei. 
The detection of broad iron emission features would argue 
against the simple, low-radiative efficiency accretion mod- 
els discussed here (Section 4.1), and require the presence 
of significant amounts of cold, reflecting material close to 
the central black holes (as in standard thin-disk accretion 
models). 
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5 PREVIOUS CONSTRAINTS ON NUCLEAR 
X-RAY EMISSION FROM M87 

5.1 ROSAT HRI observations 

Reynolds et al. (1996) and Harris et al. (1997, 1998) discuss 
ROSAT High Resolution Imager (HRI) observations of M87, 
and measure a time-averaged HRI count rate associated with 
the active nucleus of the source of ~ 0.12 count s~^. We 
have used the XSPEC software, incorporating the 1990 De- 
cember version of the HRI response matrix, to determine 
the flux of the power-law component required to produce 
the count rate observed in the 0.1 - 2.4 keV ROSAT HRI 
band. For an assumed power-law omission model with a pho- 
ton index, F = 1.4, absorbed by a Galactic column density, 
ATh = 2.5 X 10^° atom cm wc find that a normalization, 
Al ~ 1.5 X 10~^ photon kcV^^cm^'^s^^ is required to ac- 
count for the observed HRI count rate. This is in excellent 
agreement with the value measured from the ASCA spectra 
of Al = lAtol X 10"=^ photon keV-'cm-^s-\ (We note 
that accounting for the intrinsic absorption also detected in 
the ASCA spectra would imply a larger intrinsic normaliza- 
tion for the power-law component of Ai ~ 5.0 x 10~^ photon 
keV~^cm~^s~^. However, this corrected value assumes that 
the absorption is due to cold gas in a uniform screen in 
front of the nucleus, which may not be the case e.g. Allen 

6 Fabian 1997; Reynolds 1997). 



5.2 Previous analyses of the ASCA observations 

Two previous studies of M87 have also reported results based 
on the same ASCA observations discussed in this paper. 
Matsumoto et al. (1996) present results from an analysis 
using a simple two-temperature plasma plus power-law spec- 
tral model, with a fixed power- law photon index of F = 1.7, 
from which they determine a flux associated with the nucleus 
(in the 0.5-4.5 kcV band) of 1.1 ±0.2 x 10"" crgcm^^s"^ . 
This flux is ~ 2.5 times larger than the value implied by our 
best-fitting spectral model (model G). Reynolds et al. (1996) 
also present results from a study of the same ASCA obser- 
vations, in which they did not detect significant emission 
associated with the nucleus, although these authors did not 
account for the possibility of variable element abundance 
ratios in their analysis, which provides a crucial step in the 
modeling (sec also Allen et al. 1999). To further illustrate 
this, we have repeated our analysis of the power-law com- 
ponent in M87 with the element abundances linked to vary 
in the same ratio relative to their solar values (c./. spectral 
model D). The best-fit value obtained with this more sim- 
ple model, = 1888, is substantially worse than the value, 
= 1468, obtained with spectral model G. The best-fit 
parameter values (F = 1.70l";3} and Ai = l.'dtH x 10"^ 
photon keV~^cm~^s~^) are also (slightly) offset from the re- 
sults obtained with model G. Such differences demonstrate 
the need to fully account for the complex temperature struc- 
ture of the galaxy /cluster plasma and possible variations in 
individual element abundances ratios when attempting to 
constrain the power-law emission from such sources. 



5.3 Rossi X-ray Timing Explorer observations 

Reynolds et al. (1999) present further constraints on power- 
law emission from M87 using observations made with the 
Proportional Counter Array (PGA) on the Rossi X-ray Tim- 
ing Explorer (RXTE). The data presented by these au- 
thors cover the 3 — 15 keV range and constrain the 2 — 10 
keV flux of the power-law component to be < 4.1 x 10"^^ 
ergcm~^s~^ , for an assumed photon index F = 2.0. This 
upper limit is lower than the measured value of 8.7lJ;g x 
IQ-^^ ergcm-^s"^ from the ASCA data (Table 4). The 
ASCA data also require a significantly flatter photon index 
than assumed in the RXTE study. 

The best-fltting spectral model for M87 determined 
from the ASCA analysis and plotted in Fig. 2 shows that 
the extended cluster emission dominates over the power-law 
component in the ASCA spectra across virtually the entire 
energy range of the detectors. Only at energies £ > 8— 9 keV 
does the power-law component dominate the detected flux. 
The much larger field of view of the PGA collimator (~ 1 
degree^ FWHM; Jahoda et al. 1996) results in a larger frac- 
tion of the total cluster flux (--^ 5 times more) being included 
in the detected spectrum (the cluster emission extends to a 
radius of at least 4 degree; Schindler, Binggeli & Bohringer 
1999). Thus, the power-law component, as determined from 
the ASCA data, will not dominate over the cluster emission 
in the PGA spectrum below an energy of ~ 12 — 13 keV. 

Modelling the extended plasma emission in a cluster 
as cool as Virgo Cluster (Table 3) with instruments like 
the PGA, restricted to the (relatively hard) 3 — 15kcV en- 
ergy range, is difficult. The PGA spectra cannot reliably 
constrain the multiphase nature of the gas in the cluster 
core which, as this study has shown, can be crucial in con- 
straining the properties of the power-law emission. However, 
although such considerations may be relevant in interpret- 
ing the RXTE results, it remains plausible that the nuclear 
emission from M87 may simply have varied between the 
ASCA observations in 1993 June and the PGA observations 
made in 1998 January (Harris et al. 1997, 1998; Tsvetanov 
et al. 1998) 

5.4 ROSAT HRI flux-limits on other sources in 
the sample 

The other galeixies discussed in this paper have not pre- 
viously been studied in the same detail as M87 and no 
detections of point-source X-ray emission associated with 
their nuclei have been reported. However, Di Matteo et 
al. (1999a) present limits on possible nuclear X-ray emis- 
sion for the three Virgo ellipticals, based on ROSAT HRI 
imaging data. Their limits, which are defined at an en- 
ergy of IkeV, are vF^ < 6.8 x 10"" ergcm'^ s"^ (NGG 
4472), vF^ < 7.8 x 10"" ergcm'^ s"^ (NGG 4636) and 
uF,, < 1.5 X 10"" crgcm-^s"^ (NGC 4649). These lim- 
its compare to the measured ASCA fluxes, quoted in the 
same units, of uF^ ~ 6.1 x 10"" crgcm"^ s"^ (NGC 4472), 
z/F„ ~ 1.6 x 10"" ergcm-^s"^ (NGC 4636), and vF^ ~ 

3.5 X IQ-" ergcm-^s"^ (NGC 4649). 

The ASCA measurements for NGC 4472 and 4649 arc 
within the Di Matteo et al. (1999a) limits. For NGC 4636, 
however, the ASCA measurement is approximately twice the 
ROSAT limit. The Di Matteo et al. (1999a) limits are de- 
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termined by fitting an analytic King model to the observed 
X-ray surface brightness profiles and determining the maxi- 
mum additional contribution that can be made by a central 
point source. However, these limits are sensitive to complexi- 
ties in the observed surface brightness profiles and, especially 
for NGC 4636, which exhibits a complex X-ray morphology 
in its central regions, the ROSAT limits should be viewed 
with caution. 



6 THE ELEMENT ABUNDANCES IN NGC 

4636 

In contrast to the ASCA results for M87 and NGC 4696 
(Allen et al. 1999), for which the introduction of individual 
element abundances as free parameters in the fits leads to 
a more significant improvement in the statistical quality of 
the fits than the introduction of the power-law component, 
for NGC 4696, the-introduction of the power-law emission 
component provides by far the most significant improvement 
over the basic two-temperature model. Thus, in determining 
our results on the abundances of the individual elements in 
NGC 4636, we started from the two-temperature model with 
the power-law component included, which provides a of 
722.6 for 230 degrees of freedom. We then systematically 
determined the statistical improvements to the fit obtained 
by allowing the abundance of each element, in turn, to be a 
free parameter in the analysis. Having identified the element 
providing the most significant improvement, the abundance 
of that element was maintained as a free parameter, and 
the process repeated to determine the element providing the 
next most significant improvement. 

In agreement with the results for M87 and NGC 4696, 
we find that the most significant improvements in the fit 
to the NGC 4636 data are obtained by including the abun- 
dances of Si, Mg and S as free fit parameters (the measured 

value is reduced from 722.6 to 475.0, for the introduction 
of only three extra fit parameters). At this point, includ- 
ing the abundances of further elements as fit parameters 
did not lead to such significant improvements and, due to 
the already complex nature of the best-fit model, we forced 
the abundances of all other elements to vary with the same 
ratio relative to their solar values (effectively tying the abun- 
dances to that of iron, the element to which the ASCA data 
are most sensitive). We note, however, that including the 
abundances of Na and O as further free fit parameters did 
lead to formally significant improvements, with reductions in 

of ~ 20 and 30, respectively. However, due to the system- 
atic uncertainties in the NGC 4636 spectra at low energies 
(Section 3.3), which effect the measured O abundance, and 
the fact that the Na abundance fits to an un-physically high 
value (~ 5 times solar), the metallicities of these elements 
were not included as free parameters in our final analysis). 
The measured element abundances for NGC 4636, as a frac- 
tion of their solar pliotosplicric values defined by Anders 
& Grovesse (1989), arc Zfo = 0.62l[;:J-^, Zms = 0-93to ll, 
Zsi = 0.96lo;}2 and Zs = lAQ^o^l. These results are in 
reasonable agreement with those of Matsushita et al. (1997) 
and Buotc (1999). 

Scaling our measured abundance ratios to the mete- 
oric abundance scale of Anders & Grevesse (1989) we de- 
termine [Mg/Fe] ~ 0.00, [Si/Fe] ~ 0.02 and [S/FeJ^ 0.16. 



Comparing these values with the supernova yield models of 
Nagataki &: Sato (1998), our observed [Si/Fe] ratio implies 
a mass fraction of the iron enrichment due to type la su- 
pernova, MFc,sNia/MFe, total, iu the range 0.55 — 0.9 (where 
the limits cover the full range of models examined by these 
authors). For spherical type II supcrnovae, a mass fraction 
in the range .A^Fe,SNia/.^^Fe, total 

~ 0.7 — 0.85 is preferred. 
Further comparison of our [Si/Fe] constraint with the su- 
pernovae models discussed by Gibson et al. (1997) also re- 
quires A^Fc.sNia ~ 0.5 — 0.7. However, the observed [S/Fo] 
and [Mg/Fe] ratios favour a mass fraction due to type la 
supernovae of < 0.5 (Gibson et al. 1997). 



7 CONCLUSIONS 

We have presented detections of hard X-ray emission com- 
ponents in the spectra of six, nearby giant elliptical galax- 
ies observed with ASCA. The galax;ies exhibit clear, dy- 
namical evidence for supermassive (10*— a few 10^ Mq) 
black holes in their nuclei. The hard X-ray emission can 
be parameterized by power-law models with photon in- 
dices, r = 0.6 — 1.5 (mean value 1.2), and luminosities, 
Lx,i-io ~ 2.6 X 10^" - 2.1 X 10*^ ergs'S or thermal 
bremsstrahlung models with electron temperatures, kT > 10 
keV. Such properties identify these galaxies as a new class 
of accreting X-ray source, with X-ray spectra significantly 
harder than those of Scyfcrt nuclei, typical binary X-ray 
sources and low-luminosity AGN, and bolometric luminosi- 
ties comparatively dominated by their X-ray emission. We 
have argued that the hard X-ray emission is likely to be due 
to accretion onto the central, supermassive black holes, via 
low-radiative efficiency fiows, coupled with strong outfiows. 
Within such models, the hard X-ray emission originates from 
bremsstrahlung processes in the radiatively-dominant, outer 
regions of the accretion flows. (Detailed modeling and dis- 
cussion of these issues are presented by Di Matteo et al. 
1999b). 

For the case of M87, the flux of the hard component 
was shown to be in good agreement with the nuclear X- 
ray flux dotorminod from earlier ROSAT HRI observations, 
which were able to resolve knot A in the jet from the nu- 
clear emission component. We have stressed the importance 
of accounting for the complex temperature structure, intrin- 
sic absorption and variable element abundance ratios in the 
analysis of the ASCA spectra. We confirmed results showing 
that the application of such models leads to measurements of 
approximately solar emission-weighted metallicities for the 
X-ray gas in the galaxies. We also presented detailed results 
on the individual element abundances in NGC 4636. 

Future observations at high spatial resolution with the 
Chandra Observatory will be crucial in establishing the 
contributions from the various X-ray emission mechanisms 
present in elliptical galaxies and unambiguously identifying 
the origin of the hard X-ray components. Deep X-ray spec- 
troscopy with XMM and ASTRO-E will allow us to examine 
variability in the X-ray emission (which should be slower in 
sources where the X-ray emission originates from the outer 
regions of low radiative-efficiency accretion fiows than in 
typical Seyfert nuclei) and to search for broad, iron emis- 
sion features associated with the power-law components. 
The detection of such broad emission features would argue 
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against the simple, low-radiative efficiency accretion mod- 
els discussed here, and require the presence of significant 
amounts of cold, reflecting material close to the central black 
holes. 

The discovery of hard X-ray emission components in the 
spectra of nearby elliptical galaxies containing supermassive 
black holes provides important new constraints on the ac- 
cretion processes in these systems. Our results are relevant 
to understanding the demise of quasars (which could plau- 
sibly be due to a change in the dominant accretion mode in 
ellipticals over the history of the Universe) and, ultimately, 
the origin of the hard (F ~ 1.4) Cosmic X-ray Background 
{e.g. Di Matteo & Fabian 1997b). These issues, and others, 
will be explored in future papers (Di Matteo et al. 1999b; 
Di Matteo & Allen 1999). 
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